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1. Introduction 

Neumann system describes motion of a particle constrained to the n-dimensional 
sphere S n under quadratic potential. The potential is given in ambient coordinates 
q = (qi, . . . , q n +i) £ M™ +1 by the potential matrix A — diag(ai, . . . , a n+ i) as 



In the generic case where all the eigenvalues of the potential ai are different, the 
Neumann system is algebraically completely integrable and its flow can be linearized 
on the Jacobian torus of an algebraic spectral curve [U El [3] . A standard approach to 
study integrable systems is by writing down the system in the form of Lax equation, 
which describes the flow of matrices or matrix polynomials with constant eigenvalues, 
i.e. the isospectral flow 014]. Eigenvalues of the isospectral flow are the first integrals 
of the integrable system and Lax representation maps Arnold-Liouville tori into the 
real part of the isospectral manifold consisting of matrices with the same spectrum. 
A quotient of the isospectral manifold by a suitable gauge group is in turn isomorphic 
to the open subset of the Jacobian of the spectral curve [5]. 

Two different Lax equations are known for a generic Neumann system with n 
degrees of freedom: one is using (n + 1) X (n + 1) matrix polynomials of degree 2 [2] 
and the other is using 2x2 matrix polynomials of degree n j3)|5]. The (n + 1) x (n + 1) 
Lax equation was used by Audin [6j to describe the Arnold-Liouville tori for geodesic 
motion on an ellipsoid, which is equivalent to the Neumann system. 

In contrast to the generic Neumann system, the special confluent case in which 
some eigenvalues of the potential coincide has not received much attention. In this 
paper we will consider the confluent case with two of the eigenvalues being the same 
and the Neumann system having an additional S 1 rotational symmetry. We will 
show that the confluent Neumann system is algebraically completely integrable and 
that its flow can be linearized on the generalized Jacobian of a singular algebraic 
curve. We will describe the symplectic reduction of the S* 1 action, which will yield 
an alternative description of the Rosochatius system as a symplectic quotient of the 
confluent Neumann system. Mechanically speaking, Rosochatius system can be seen 
as a Neumann system on a rotating sphere. Inspired by this result we will describe 
general Rosochatius system as a reduction of confluent Neumann system with all 
eigenvalues of the potential being double. This result combined with the proof of 
integrability of the confluent Neumann system will also give an alternative proof of 
the algebraic integrability of the Rosochatius system. 

When applying the 2x2 Lax equation to the confluent case, the resulting Lax 
equation in fact describes the Rosochatius system and not the confluent Neumann 
system [71IH1 E] ■ We will therefore use (n+1) x (n+1) Lax equation, where the resulting 
spectral curve is singular in the confluent case. Following the standard procedure and 
normalizing the spectral curve results in the loss of one degree of freedom. In order 
to avoid that, we will use the generalized Jacobian of the singular spectral curve to 
linearize the flow as in [TO] ITT] . The generalized Jacobian is an extension of the 
"ordinary" Jacobian by a commutative algebraic group (see [U] for more detailed 
description). In our case the extension will be the group C* that corresponds to the 
rotational symmetry of the initial system. Generalized Jacobian was used by others 
to linearize the flow of other integrable systems with rotational symmetry for example 




i=l 
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spherical pendulum or Lagrange top O HH HU [10] . Our study however will give a 
more detailed description of the relation between symplectic reduction and algebraic 
reduction from generalized to ordinary Jacobian. 

After the introduction, Hamiltonian reduction and Liouville integrability of 
the confluent Neumann system are discussed in section [2l In section [3] we study 
(n + 1) x (n + 1) Lax equation and corresponding isospectral manifolds. Our main 
result is formulated in theorem 13.51 and describes the relation of Arnold-Liouville tori 
to the generalized Jacobian of the singular spectral curve. As a corollary the complete 
algebraic integrability and Liouville integrability of the confluent Neumann system will 
follow. We conclude with proposition I3.8[ which describes the bifurcation diagram of 
the energy momentum map in terms of algebraic data. 

2. Hamiltonian description 

The Neumann system describes a particle on a sphere (of arbitrary dimension) under 
the influence of quadratic potential. We can write it as a Hamiltonian system on 
the cotangent bundle of the sphere T*S n with canonical symplectic form lu c and the 
Hamiltonian H given in ambient coordinates (q,p) 6 R n+1 x R™ +1 as 

H(q,p) = ±(\\p\\ 2 + (q,Aq}). 

The potential is given by a positive definite linear operator A on R" +1 . For simplicity 
we will assume that A is diagonal with positive eigenvalues a^, A consequence of 
positivity is that the Hamiltonian is proper and the energy level sets - and hence 
Arnold-Lioville tori - are compact. The equations of motion in Hamiltonian form are 

q=P (1) 
p = - Aq + eq 

where e = ||p|| 2 + c is chosen so that \\q\\ = 1 and the particle stays on the sphere. 
2.1. Reduction of the symmetry 

Let us consider the confluent case where all the eigenvalues a,i of the potential A are 
distinct, except a n = a n +i- The action ip' of S 1 — SO(2,M) on S n , given by rotations 
in the q n , q n +i plane, leave the potential invariant and can be lifted to the symplectic 
action 

p : T*S n xS 1 ^ T*S" 

on the cotangent bundle that leaves the Hamiltonian H invariant. We would like to 
reduce this action and describe the resulting reduced system in more detail. 

Let K be the moment map for the lifted action ip. The map K is a real function 
on T*S n as so(2)* ~ R and is the angular momentum for the rotations in q n ,q n +i 
plane 

K{q,p) = q n p n +i - q n +iPn- 

Since the action ip is not free, the description of reduced system is more complicated. 
Let us denote with F the set of points on the sphere, where the action <p' is not free. 
This is precisely the fixed point set given by the condition 
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The set F is a codimension-2 great sphere on S n . Note that on F, the value k of the 
moment map K equals 0. The set of regular points S^ eg = S n — F is an open subset 
of S n on which the action is free. 

The reduced system is defined on symplectic quotients (M k ,u> k ), which are 
parametrized by the value k of the moment map K. We will use the operator j j k to 
denote the symplectic quotient. By definition 

M k = T*S n // k S l := K-\k)/S l 

and uj k is defined with Tr*uj k = Uc\K- 1 (k)> where 7r is the quotient projection. For 
k ^ the fiber K~ 1 (k) is a corank-1 sub-bundle of the cotangent bundle over the set 
of regular points T*S™ eg . Since the action is free on K^ 1 (k), the standard result for 
lifted actions gives 

where uj p is the magnetic term coming from the curvature of the mechanical 
connection. In our case the mechanical connection is flat and oj p = [15] . The 
quotient manifold S^g/S 1 is a dimension-(n — 1) open half sphere S"^ 1 as we can 
see if we introduce cylindrical coordinates (q±, . . . , q n -i, r (cos ip, sin </?)). The set S™ eg 
is given by the condition r / and the quotient S^g/S 1 can be parametrized by 
(qi, . . . , q n -i, r) G S n ^ 1 for r > 0. We can also see directly by using cylindrical 
coordinates that the perturbation uj p of the canonical symplectic form is zero. 

For k = the description of Mq is more complicated, since the action on if~ 1 (0) 
is not free, and Mq is not a manifold. We will remedy this by considering a singular 
double cover of the reduced phase space instead. 

The set if -1 (0) is not a sub-bundle of T*S n , because at fixed points r = the 
fiber has full rank and over 5" eg the fiber has codimension 1. However, we can still 
pass the quotient on the base manifold because the action is a cotangent lift. The 
quotient space S n /S 1 is a dimension- (n — 1) closed half sphere (S™~ 1 ) c and can be 
parametrized by cylindrical coordinates (qi, . . . , q n -i, r) S S n , with r > 0. The 
symplectic quotient Mo is a cotangent bundle over S n /S 1 , in the sense that the fiber 
at a singular point with r = is a closed half plane of "positive differentials" (the 
differential, which is positive at the directions normal to the boundary of S n /S 1 ). In 
that sense we can write 

(M ,c ) = (T*(5r 1 ) c ^ c ). 

To conclude the construction of the reduced system we have to calculate the 
reduced Hamiltonian. This can be easily done using the parametrization with 
cylindrical coordinates (q,p) = ((qi,.. . ,g„_i,r), (pi,... ,p n -i,p r )) 

1 

H(q,p) = - I 22 (p^ + a t q^) + pi + p 2 n+1 + a n (ql + g£ +1 ) 
\t=l 

1 / n_1 k 2 
= 2 ^ p2% + ^ + ^2 + aiq2 i + a " r2 

\i=l i=l 

= \ (||P|| 2 + V k (q) + {q,Aq)) = H r (q,p), 

where A = diag(ai, . . . , a n ). The additional term V k (q) = ^ corresponds to the 
centrifugal system force, induced by the rotation. 
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To simplify the description of the reduced system, we will use the whole sphere 
instead of only half of it. The half sphere (S'" _1 ) c can be also viewed as a 
quotient of S 1 ™ -1 by the group Z 2 acting with reflections r t— > — r. The reduced 
Hamiltonian H r (q,p) can be lifted on (T*S n , u c ), since it only depends on r 2 and 
pi and is invariant for Z2 action. We have thus constructed a Hamiltonian system 
(T*S n ^ 1 , u c , H r ) such that its reduction coincides with the reduced Neumann system. 

Theorem 2.1. The singular symplectic quotient (Mk,Uk,H r ) of the Neumann 
Hamiltonian system (T* S"^ 1 ,u c , Ha) with potential given by the matrix A = 
diag(ai, . . . , a n , a n ) is isomorphic to the quotient by the 1% action of the perturbed 
Neumann system (T*S n , u> c , H^ + V&) with potential matrix A = diag(ai, . . . , a n ). 

The symplectic quotient of confluent Neumann system is in fact a special case 
of Rosochatius system, which was studied in [Tj for example. A general Rosochatius 
system is a Hamiltonian system on T*S n with potential 

We can generalize the procedure from the above to the case where potential matrix 
has all the eigenvalues double A = diag(ai, ai, a^, 0,2, ■ ■ ■ , a n +x, a n +i)- This gives a 
mechanical interpretation of Rosochatius system as Neumann system on a rotating 
sphere. 

Corollary 2.2. The symplectic quotient of the confluent Neumann system with 
potential matrix A — diag(ai, 01, 02, 02, • • ■ , ctn+l> o>n+i) by the group of rotations 
(S 11 )" - " 1 " 1 is isomorphic to the quotient by the group of reflections Z 2 l+1 of the 
Rosochatius system on T*S n given by the potential The coefficients fcj of the 

rational part are the values of the angular momentum of rotations in corresponding 
eigenplanes of the potential matrix A. 



2.2. Integrability 

Hamiltonian system with Hamiltonian H on the symplectic manifold of dimension 
2n is called integrable if there exist n functionally independent pairwise Poisson 
commutative first integrals, one of which is H. To show that the system is completely 
algebraically integrable, we will prove that the level sets of the first integrals are real 
parts of the extensions of Abelian varieties by C* (see [3] for a definition of complete 
algebraic integrability) . 

The integrability of the confluent Neumann system is a consequence of the 
integrability of the generic Neumann system, which is well knownpQ Hfi]. The first 
integrals for the confluent case can be obtained by taking the limit a n — > a„+i on 
Uhlenbeck's integrals for generic Neumann system: 

Ff = q? + E {qiP ° = qiPl? ■ (3) 

... CLj CLi 

3+1- 

A set of commuting integrals, equivalent to JH}, tends to a set of commuting integrals 
for the confluent case 

{Ff, F 3 n _ x ,F° n + Ff l+1 , \{a n - a n+1 ){F^ - F 9 n+l )} - 

^{F 1 ,...,F n _ 1 ,F n ,K 2 } 
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when taking the limit a n — ► a n+ i. The set of commuting integrals for the confluent 
Neumann system is given by 

" " "< 

jp 2,2 , \ - {qjPn - qnPjf + (qjPn+1 ~ Qn+lPj) 2 . , 

= 9„ + «n+l + 2^ ~7T~^~a ^ > 

ilj tl n 

J<11 J 

K 2 = {q n p n +i - q n +lPn) 2 - 

Note that K is angular momentum for the rotations in q n ,qn+i plane. We can also 
verify that the integrals in |f4]) are not independent but satisfy the same relations as 
Uhlenbeck's integrals in the generic case 

n 
i=l 

and that the Hamiltonian H can be expressed as a linear combination of Fi and K 2 

n 

atFi + K 2 = 2H. 

i=l 

The Poisson brackets of Fi are continuous functions of at so the commutativity of 
the integrals is preserved when taking the limit a n — > a„+i. Commutativity of K 
with Fi also follows from the fact that F{ are invariant for S 1 action, generated by 
K. To conclude the proof of integrability one needs to verify the independence of the 
integrals Q (up to relation = !)• The commutativity of the integrals and K 2 
also follows from the AKS theorem, if we write Fi an K 2 as invariant functions of Lax 
matrix on appropriate loop algebra. This is a standard way to prove commutativity 
of Uhlenbeck's integrals for the generic Neumann case 
Let us combine all the first integrals in a map 

F EM : T*S n -» E n 

(q,p) ^ (F 1 ,...,F n _ 1 ,K) 

we will call energy momentum map. The fundamental property of integrable systems 
is that the level sets of energy momentum map Fem are n-dimensional tori, on which 
the flow can be linearized. We will also use the complexified version of Fem, defined 
on {T*S n ) c . 

Theorem 2.3. The confluent Neumann system (T* S n ,ui c ,H) is algebraically 
completely integrable system. 

We have seen that by taking the limit a n — > a n+ i, {F\, . . . , F n _x, K 2 } is a set 
of n commuting first integrals for symmetric Neumann system and we will show later 
in 13.31 that they are functionally independent. We will also prove the part about 
algebraically complete integrability in subsection 13.31 

From the fact that the symmetric Neumann system is integrable also follows that 
its symplectic quotient is integrable. This gives alternative proof of the integrability 
of Rosochatius system, which is a symplectic quotient of the symmetric Neumann 
system. 

Corollary 2.4. The Rosochatius system (T*S n , u c , H + V r ) with V r = £ ^ is 
algebraically completely integrable system. 
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3. Lax representation of the confluent Neumann system 

In this section we will use Lax equation for (n + 1) x (n + 1) matrix polynomials to 
study confluent Neumann system with n degrees of freedom. We will show that the 
flow of the system can be linearized on the Jacobian of the singular spectral curve 
and that the system is completely algebraically integrable. This will also conclude 
the proof of Liouville integrability and yield the description of bifurcation diagram for 
energy momentum map. 

3.1. Lax equation in Ql(n + 1,C) 

Let us write the confluent Neumann system ([I]) as an isospectral flow of matrix 
polynomials in loop algebra gl(n + 1, C). We will use the Lax equation introduced 
by Moser [I]. The loop algebra gl(n + 1, C) consists of Laurent polynomials with 
coefficients in gl(n + 1, C), and can be written as a tensor product gl(n + 1,C) <g> 
C[A, A -1 ]. Consider complexified Neumann system (TJ on a subspace (T*S n ) C 
C" +1 x C n+1 defined by constraints J2lf = 1 an d Y^QiPi = 0- We introduce Lax 
matrix polynomial from gl(n + 1, C) 

L{X)=AX 2 +qApX-q(g>q; q,peC n+1 (5) 

where A is the potential matrix of the Neumann system and (q,p) E (T*S n ) c . 
Neumann system either generic or confluent can be written as Lax equation 

±L(X) = [M(X),L(X)]. (6) 

with M(A) = AX + qAp = (A£(A))+ (the subscript + denotes the polinomial part of an 
element of gl(n + 1, C))[2]. The flow of ((H) is isospectral as it conserves the spectrum 
of the matrix L(X). This means that the characteristic polynomial P(A,/i) and the 
corresponding afflne spectral curve, defined by equation 

P(A,/i) = det(L(A) -m) = 

do not change along the flow and can be expressed only with the values of the first 
integrals of Neumann system. In order to avoid unnecessary singularities at the 
infinity, the affine curve is completed in the total space of the line bundle Opi (2) over 
P 1 , which is given by the transition function (A, fi) (A -1 , A~ 2 u) . The completion of 
the affine spectral curve in Opi(2) is called the spectral curve of L(X) and denoted by 
C m . Note that for any given L(X) there are also naturally defined a map A : C m — > P 1 
and a section [i of the bundle A*0pi (2) apart from the spectral curve C m . 

The spectral curve of L(X) is hyperelliptic and we can see that by introducing 
new variables x = A~ 2 /i and y = AJ}^ (dj — x), where a, are the eigenvalues of the 
potential matrix A. The equation we obtain is: 

n+1 

y 2 = Q(x)H(ai-x), (7) 

i=i 

where Q(x) is a polynomial of degree n [H [6]. The coefficients of Q(x) are first 
integrals of the confluent Neumann system and if we write Q(x) in terms of Lagrange 
interpolating polynomials over the points ai, . . . ,a n ,a n using notation with partial 
fractions 

Q(*) = Ui<h-*)[-^ + T J ^ + 'L — )> W 
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the coefficients we obtain are the integrals ((!]) we met before . 

In the confluent case, where a n — a n +i, the product Yl i a i ~ x ) has a quadratic 
factor (a n — x) 2 and the spectral curve has a singular point S given by (x, y) — (a n , 0) 
or (A, /i) = (oo,a n ). The singularity S is a double point for K ^ and a cusp for 
K = 0. The smoothness of the spectral curve is closely related to the regularity of the 
matrix L(X). Recall that a matrix B £ gl(r, C) is called regular if all the eigenspaces 
of B are one dimensional. 

Proposition 3.1. Let C be the spectral curve of matrix polynomial L(X) and a £ P 1 . 
If all the points A (a) £ C are smooth, then the matrix L(a) is regular. 

For proof see [5j [10]. The value of L(A)A~ 2 at A = oo is the matrix A and the 
singularity at S is a consequence of the fact that A is not regular when a n — a n+ \ . 
Let C be normalization of the spectral curve C m , which is described by 

n-l 

w 2 = Q(x) J[(ai-x), (9) 

i=l 

where w = y/(a n — x). We will call C the normalized spectral curve of L(X). There 
is a map tt : C — > C m that is biholomorphic everywhere except at the inverse image of 
the singular point S. The inverse image tt~ 1 (S) consists of two points {P + ,P-} for 
and a point Po for K = 0. In case when K ^ 0, the curve C m is obtained from 
C by identifying the points {-P+, P— } into the singular point S. The singular curve C m 
can be described as a singularization of C given by modulus m (see [12] for details) . 
The modulus is m = P+ + P_ for K ^ and m = 2F for if = 0. 
Remark 3.2. Note that we will only resolve the singularity at S, which is "generic" 
in the sense that it appears for all the values of the energy-momentum map. So the 
curve C can still be singular for some values of the energy- momentum map. 

Finally we find the genus of C and C rn from the fact that the curves are 
hyperelliptic and from the degree of the polynomials in ([7]) and (J9j) . We obtain 
g{C) = n — 1 for normalized spectral curve C and the arithmetic genus g a {C m ) = n for 
the singular curve C m . 

3.2. Isospectral manifold of matrix polynomials 

We have seen that the Neumann system satisfies Lax equation and that the spectral 
curve depends only on the first integrals and that all the first integrals are encoded 
in the spectral curve. The level sets of (complexified) energy momentum map Fem 
lie in the set of matrices L(A) with fixed spectral curve C m . It is therefore essential to 
describe the set of matrix polynomials with a given spectral curve. 

Let P(A, jj) be a spectral polynomial for some matrix polynomial L(X) £ fll(r, C)). 
Denote with M. p the subset of all the elements of gl(r, C) with the same characteristic 
polynomial P (this also fixes the spectral curve C) 

M P = {L(X); det(L(A) -/i) =P(A,/x)}. 

All L(X) £ A4p have the same spectral curve C, which is defined by P(X,/i) = 0. 

While the characteristic polynomial and thus the eigenvalues of L(X) are fixed by 
the flow, the eigenvectors and eigenspaces change. Let us define a map 

such that £l(a)((A, A*)) is one dimensional eigenspace of the matrix L(X) with respect 
to the eigenvalue [i. If the spectral curve is smooth, than by proposition 13.11 all the 
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eigenspaces of L(A) for any A are onedimensional and the map £ is well defined. The 
map £l(a) defines a line bundle on C and its dual is called eigenvector line bundle or 
shorter eigenline bundle. We will denote eigenline bundle by Cu\y By construction, 
the eigenline bundle is a subbundle of the trivial bundle CxC. One can see by using 
Riemann-Roch-Grothendick theorem that the degree (Chern class) d of the eigenline 
bundle Cl(x) equals g + r — 1 where g is the genus of the spectral curve C. 

The only condition for £z,(A) to be defined is that the eigenspaces of L(X) are one 
dimensional for all but finite number of points on C. If the spectral curve is singular, 
than the map £l(A) can be defined on the set of points C rn — N C C m where the 
matrix L(X) has one dimensional eigenspace. If the set N is finite, the map £l(a) can 
be extended as a holomorphic map and eigenline bundle Lux) can be defined on the 
normalization C of the spectral curve. Note that by the proposition 13.11 the set N is a 
subset of singular locus of the spectral curve C m . 

One can define the eigenbundle map 

e:M P -> Pic d (C) 
L(X) •-> [£ LW ], 

from Mp to the Picard group Pic(C) of isomorphism classes of line bundles on the 
normalized spectral curve C. The subset Pic d (C) consists of classes of line bundles with 
given Chern class d. The set Pic d (C) is a copy of the zero degree Picard subgroup 
Pic .(C), which is in turn isomorphic via Abel-Jacobi map to the Jacobian Jac(C) 
of C. The map e assigns to each matrix polynomial L(X) the isomorphism class of 
its eigenline bundle and thus encodes the flow of Lax equation. The map e is not 
surjective since eigenline bundles cannot lie in the special divisor 6 on the Jacobian. 
The map e is neither injective since the class of eigenline bundle defines the matrix 
polynomial only up to conjugation by the gauge group ¥Gl(r, C). The space we have 
to consider is the quotient space Mp/PGl(r,C) and it was shown in [5] that if the 
spectral curve C is smooth, the space Mp/PGl(r, C) is isomorphic as an algebraic 
manifold to the Zariski open subset Jac(C) — 9 of the Jacobian of the spectral curve 
C. The isomorphism is given by the eigenbundle map e. 

As the leading coefficient in L(X) is preserved by the flow, we will use the 
closed subset of Alp of matrix polinomials with fixed leading coefficient. Let 
L(X) = AX 1 + A^iX 1 - 1 + ...+A and let A be fixed. We denote 

Mp- = {L(X) G M P ; lim L(X)/X l = A}. 

A — >oo 

The action of the gauge group PGZ(r, C) on Mp reduces to the action of the stabilizer 
subgroup FGa < VGl(r,C) of A. The quotient Mp by PGa is again isomorphic to 
the Jacobian 

Mj,/fG A ~ Jac(C) - 9. 

Remark 3.3. If the spectral curve C is smooth at infinity, the matrix A has to be 
regular by proposition 13. 11 As a consequence, the stabilizer group of A is the product 
VGa = (C*) s_1 x C _s where s is the number of distinct eigenvalues of A. If A is 
not regular, the dimension of the stabilizer group PGa is larger. In our case, when 
A = diag(ai, . . . , a n , a n ) the group PGa is the product 

(C*)"- 2 x Gl(2,£). 

Let us assume for one moment that the matrix A is regular. We have seen before 
that the level set Mp is an extension of Jac(C) — 6 by an Abelian algebraic group 
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VGa- As it happens, the generalized Jacobian is also defined as an extension by 
Abelian algebraic group and it was shown in [10] that 

Mp ~ Jac(C m ) - 

for a suitable choice of modulus m. The chosen modulus is the effective divisor 
consisting of infinite points Pi on the spectral curve 

where the coefficients m, are multiplicities of the eigenvalues /u(Pj) of A — L(oo). 
The isomorphism A4p — ► Jac(C m ) — 8 is given by the eigenbundle map e m to the 
generalized Jacobian. Let S — \m\ be the set of points in C that are mapped to the 
singular point in C m . The line bundle on the singular curve C m is given by a divisor 
on C that does not intersect the singular set S, To give a line bundle on C m is thus 
enough to give a section of a line bundle on C that has no zeros or poles in S. We can 
then extend the map e : Mp — > Pic d (C) to a map 

e m : M P -> Pic d (C m ) 

by choosing a section of e(L(A)) uniformly on Alp that does not have any zeros or 
poles in the infinite point. This can be done by appropriate normalization. It was 
shown in |10[ fTT] that e m gives an isomorphism from the isospectral space A4p to the 
Zariski open subset Jac{C m ) — of the generalized Jacobian of the singular spectral 
curve, given by modulus m. 

Remark 3.4. The above results give orientation about the expected number of degrees 
of freedom of isospectral flows. We see that the upper limit is the arithmetic genus of 
the singularization of the spectral curve, which depends on the degree I and the rank 
rofL(A). 

In the case of confluent Neumann system, the matrix A is not regular and 
corresponding spectral curve is singular. We will show later that the Lax flow also 
preserves part of the "derivative" of L{\) at the singular point. We will restrict the 
space Aip further by fixing a specific block of the lower term Ai-\. The resulting 
isospectral set will again be isomorphic to the open subset of the generalized Jacobian 
[TT] . Note that we have to assume regularity of the fixed block of the lower term Ai-i 
in order to have the isomorphism. 

3.3. Proof of the integrability 

In order to prove the Liouville integrability of the confluent Neumann system, we have 
to show that its first integrals {Fi, . . . , F n —i,K} Poisson commute and that they are 
functionally independent. We already know that the integrals commute from section 
[2j We also know that the level sets of the first integrals - the Arnold-Liouville tori - lie 
in the isospectral manifold A4p. We will use the eigenbundle map e m to map Arnold- 
Liouville tori to the real part of the generalized Jacobian of the singular spectral curve 
C m and show that this map is a (Z 2 )™ -1 covering. This will prove complete algebraic 
integrability of confluent Neumann system and independence of the first integrals will 
follow. 

Let E n be the eigenspace of the double eigenvalue a n , which is spanned by the 
unit vectors e n and e n +\ ■ The behavior of L(X)\e„ near infinity is given by the n, n + 1 
block F of the matrix q Ap, which is in fact conserved by the isospectral flow. The 
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block F depends only on the angular momentum K — q n p n+ i — q n+ ip n and is a regular 
matrix 

K 



F = 



K 



of rank 2 if K ^ 0. The isospectral flow given by J6J conserves both matrix A and 
the block F, therefore it is reasonable to consider isospectral manifold of all matrix 
polynomials with this data fixed 

Mp' F = {L(X)eM P ; L(oo)=A, w En ° i'(w)k = F } 

where pr^„is a projection to E n , and the values of L(X) and L'(X) at infinity are defined 
by the limits L(oo) := lim A ^ 00 (A" 2 i(A)) and L'(oo) := lim A ^oo d(A~ 2 L(A))/d(A _1 ). 
On Alp' acts the subgroup PGa,f < PGa of matrices that stabilize F as well. We 
can write 

PGa.f (C*)"" 1 x G F 

where Gp is the stabilizer subgroup of F in GZ(2, C). The group Gf — C* consists of 
matrices 

cos ip sin 
— sin tp cos <p 

where r s (0, oo) and </? S [0, 27r). We will describe the isospectral manifold M.^: F as 
an open subset of the generalized Jacobian of singularized spectral curve. 

Theorem 3.5. Let / G C" be the value of Fem such that the normalized spectral 
curve C is smooth and K ^ 0. Let denote by Ta — (C*)™ _1 the subgroup o/PGa of 
diagonal matrices G — [gi.j] with g n>n = g n +x >n ±i, 

(i) The complex level set FemU) * s a covering of the isospectral manifold M P ' /Ta- 
The fiber of the covering is the same as the orbit of the group (Z2)™ -1 generated 
by reflections ^ —qi for 1 < i < n — 1. 

(ii) The isospectral manifold M P V ' F /Ta is isomorphic to the open subset of the 
generalized Jacobian of the singular spectral curve C m ' given as a singularization 
of smooth spectral curve by modulus m = (00, iK) + (00, —iK). 

(Hi) The flow of K generates the fiber C* of the extension C* — > Jac(C rn ) — ► Jac(C). 

(iv) The flow of H on the generalized Jacobian Jac(C rn ) is linear. 

Proof. If K ^ then F is a regular matrix with 2 different eigenvalues ±iK. By the 

.A F 

result in [IT] the isospectral manifold M P ' is isomorphic as an algebraic manifold to 
the Zariski open subset J(C m ) — 9 of the generalized Jacobian of the singularized 
spectral curve, given by spectral polynomial P and modulus m' = A (00) = 

same theorem asserts that the following 

diagram commutes 

FG a ,f Mp F Mp F /¥G A>F 



*- (C*) n ^ Jac(C m <) - 6 > Jac(C) - 6 
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If we singularize C only by modulus m = (oo,iK) + (oo, — LK"), we can insert 
e m : Alp' F /(C*)" _1 — > Jac(C m ) — 9 into above diagram 



M 



A.F 



G, 



■M A ' F /T A 



■ Mj,> F /m A>F 



Jac{C m i) — 9 Jac(C m ) — 9 s~ Jac(C) — 9 



C* 

We only have to prove that the fiber F^ M (f) consisting of matrices of the form ([5]) 
is a covering of the quotient M^ F /T A . By using the Lax matrix J5]) we gave a 
parametrization 

J A : (q, p) h-> L(X) = AX 2 + q Ap X - q® q 

of the quotient Mp' F /T A by (q,p) G F^uif) C (T*S) C . First note that the map 
J A : (T*S n ) c — * A^p ,i? is an immersion. We will show that any orbit of intersects 
the image of J A only in finite number of points. To explain how FG a ,f acts on the 
Lax matrix ([5]), note that an element g G PGa,f acts on a tensor product x (8 y of 

!,)/£ C" +1 

g:x®y^ (gx) <g> ((g _1 ) T y) 

by multiplying the first factor with g and the second with (g _1 ) T . The subgroup of 
FGa,f, for which the generic orbit lies in the image of J A , is given by orthogonal 
matrices 

0(n, C) n PG a ,f (Z 2 )"" 2 x G F n 0(2, C). 

There are special points in the image of J A that have a large isotropy group (take 
for example qi — #y and Pi = Sik, k ^ j, where the isotropy is (C*)" -2 ). But 
the intersection of any orbit with the image of J A coincides with the orbit of 
(Z 2 )™~ 2 xG f n 0(2, C). If we take the torus T A < FG a ,f consisting of diagonal 
matrices with g n . n = g n +i,n+i> so that T^flGj? = {Id}, the orbits of will intersect 
image of J A only in the orbit of the finite subgroup (Z2)™~ 2 . We have proved that 
the level set of Lax matrices L(X) — J A (q,p) with fixed characteristic polynomial P is 
an immersed submanifold in A4p' F that intersects the orbits of torus ~ (C*)™ _1 
in only finite number of points and is thus a covering of the quotient M. p ' F /Ta- 

The group Gf is the complexification of the group of rotations in q n , q n +i plane 
and is generated by the Hamiltonian vector field of K. This proves the assertion fm|) . 

To prove the assertion (|Tv|) . note that the matrix polynomial M(X) in Lax equation 
iJH) is given as a polynomial part R(X, L(X))+ for a polynomial R(z, w) — zw. It is well 
known that such isospectral flows are mapped by e m to linear flows on the Jacobian 
Jac(C m ) (see [4] for reference). □ 
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Taking into account the real structure on C m , the Arnold-Liouville tori can be 
described as a real part of the generalized Jacobian. 

Theorem 3.6. For K ^ and C smooth, the Arnold-Liouville tori are (Z^)" -2 
coverings of the real part of the generalized Jacobian. The rotations generated by K 
are precisely the rotations of the fiber S 1 in the fibration S 1 — > Jac(C m ) R — > Jac(C) R , 
which is the real part of the fibration C* — > Jac(C m ) — * Jac(C). 

Above theorem gives us an algebraic way to describe symplectic quotient 
T*S n //kS , Algebraically T*S n is a covering of the relative generalized Jacobian, 
which is a disjoint union 

Jac R = U Cm Jac R (C m ) 

over the space of curves C m corresponding to the real values of energy momentum 
map. The symplectic quotient of Jac^f/kS 1 is then the relative Jacobian over the 
space of normalized spectral curves with fixed value of K 

U Cm ;K=fc^ac R (C). 

Corollary 3.7. The complex level set of (Fx, . . . , of the symplectic quotient of 

the confluent Neumann system is a (Z 2 )™~ 2 covering of the quotient M^,' F /PG a,f ■ 
The manifold A4 P ' F /¥Ga,f is isomorphic to the open subset Jac(C) — of the 
Jacobian of the normalized spectral curve. 

This result agrees perfectly with the results obtained previously for the 
Rosochatius system [8J. 

Proof of the theorem HOI Note that the eigenvalues of F are ±iK. Note also that the 
value of fx at the points P± equals to the eigenvalues of F, so P± = (oo, ±\K). On 
C there is a natural real structure J induced by the conjugation on (A, fi) £ C 2 . The 
points P± that are glued in the singular point form a conjugate pair P± = JP T . If 
we follow the argument in [l4j we can find the real structure of the fiber C* in the 
extension C* — > Jac(C m ) — > Jac(C). Note that Pic(C) is defined as the space of all 
divisors modulo divisors of meromorphic functions on C, whereas Pic(C m ) is given by 
the divisors on C that avoid P± modulo meromorphic functions on C m . So the fiber C* 
is given by meromorphic functions on C modulo meromorphic functions on C m . Since 
we obtained C m by gluing two points P± , a function / on C defines a function on C m 
if f(P+) = f(P-) or equivalently /(P+)//(P_) = 1. For a divisor of any function / 
on C, the number z = /(P+)//(P_) S C* determines its class in the Picard group 
Pic(C m ). So if P± — JPp, then the real structure on the fiber C* is given by the map 

f{p+) JjlPT) _ Wl = i 

f(P-) ^ f(JP-) f(P+) z 

and the real part of C* is the unit circle S* 1 given by zz = 1. In contrast, when the 
singular points are real P± = JP±, the real structure on C* is given by the conjugation 
and the real part of C* is R*. □ 
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3.4- Bifurcation diagram 

We can use the normalized spectral curve C to describe the singular locus of energy 
momentum map. 

Proposition 3.8. The vector (/i, . . . , f n -i, k) £ R" is a regular value of the real 
momentum map 



if and only if the normalized spectral curve C is smooth. The singular locus of the map 
Fem consist of 

• hyperplanes Fi = 

• zero level set of the discriminant of Q{x) from ^ 

• the codimension 2 hyperplane defined by K = and H = |aj 

Proof. The if part follows directly from theorem 13.51 If C is smooth and K ^ 0, then 
the level set of Fem is locally isomorphic to the real part of Jac(C m ) — 6, which is 
of dimension n. Hence the rank of the differential of Fem is also n. If K = and 
C is smooth, the level set of (F%, . . . , F„_i) is locally isomorphic to the real part of 
the isospectral manifold Mp' F /FGa,f, which is in turn isomorphic to the real part 
of Jac(C) — 0. Since the dimension of Jac(C) is n — 1, the rank of the differential of 
[F\, . . . , -F n _i) is also n — 1. The integrals Fi are invariant to the rotations generated 
by K and therefore their Hamiltonian vector fields Xp x , . . . ,Xp n _ 1 are independent 
from Xk- We are left to show that if Xk — the normalized curve C is not smooth. 
It is easy to see that the case Xk = appears only if Xh — but then the rank of 
Xf x , ■ ■ ■ , Xp n _ 1 is not full and the curve C has to be singular. 

To prove the only if part let us consider case by case the components of the 
singular locus. The curve C is singular if and only if the polynomial Y[( a i ~ x)Q(x) 
has a double root. This can happen in two cases 

(i) di is a zero of Q(x), this happens when Fi = 

(ii) Q(x) has double zero, this happens if the discriminant of Q is zero. 

The hyperplanes Fi = are singular, since for the points (<?,p) with q L = Pi = the 
differential di 7 !; = 0. The proof that the discriminant of Q is singular can be found 
in [6j and I will omit it here, because it is very specific and beyond the scope of this 
article. □ 

Remark 3.9. For values of a n > dj for some 1 < j < n the bifurcation diagram has a 
singular "thread" of focus-focus singularities defined by values K = and H — \a,j. 
This would suggest the presence of nontrivial monodromy. Indeed for two degrees of 
freedom, the singular level set corresponding to the isolated singular value is a union 
of two spheres with two pairs of points identified. By the general result in [17] it 
follows that the monodromy is nontrivial and equals 



Fem — {Fi , . . . , F n -i, K) 




with e = 2 being the number of spheres in the singular level set. This can also be 
checked by direct calculation. 
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K 



(a) 01 > (J2 



(b) 0,2 > a\ 



Figure 1: Image of the energy- momentum map (A", 2H) for confluent Neumann system 
with potential matrix A — diag(ai, a 2 , 02). Regular values lie in the shadowed area 
while the solid curves contain singular values. The points (0, ai) and (0,a2) are the 
images of fixed points. 

The image of the energy momentum map (K, 2H) for confluent Neumann system 
with two degrees of freedom is depicted in figure [TJ Regular values lie in the shadowed 
area while the solid curves contain singular values. The points (0, a±) and (0,02) are 
the images of fixed points. Note that when passing from 0,2 > ai to a\ > a 2 , the 
pair of lines becomes imaginary and only their intersection - isolated focus-focus point 
(0, a\) - remains real. 

3.5. Note on the case K = 

We have seen in the previous section that the case K — is significantly different 
from the generic case K ^ 0. The problem lies in the following observation. The 
space of hyperelliptic curves that appear in the description of Neumann system is 
parametrized by K 2 and not by K. As a consequence, the relative generalized Jacobian 
is degenerated for K = 0. The phase space of complexified Neumann system is 
therefore "folded" into relative Jacobian by the map K — » K 2 . The map between 
original phase space and its image in relative Jacobian is singular at K = 0. It is 
therefore illusory to expect that we can describe the whole phase space including the 
fiber K = by algebro-geometric methods. Different approach has to be considered 
that would study the "fold" given by K 2 in more detail. This is to be covered in our 
future work. 

4. Conclusions and discussion 

We have proved the algebraic integrability of the confluent Neumann system by 
proving the theorem 13.51 which describes Arnold-Liouville tori in terms of the 
generalized Jacobians of singular spectral curve. We performed the reduction of the 
rotational symmetry and established a firm relationship between symplectic reduction 
and desingularization of the spectral curve (corollary I3.7|l . Most of our results 
very likely generalize to any Moser system arising from the rank 2 perturbations 
of a fixed matrix with a double eigenvalue. From our work and previous examples 
1 141 [T3] it appears that there generally is a relation between the rotational symmetry 
and singularities of spectral curves. We have exposed this relationship explicitly 
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in our case and have seen that the reduction of S 1 symmetry reveals itself in 
the algebraic description as a reduction from generalized Jacobian of the singular 
spectral curve to the "ordinary" Jacobian of the normalized spectral curve. One can 
say that the desingularization of the spectral curve corresponds to the symplectic 
quotient. Unfortunately, this relation is not a general phenomenon as we can see 
when considering the case K = 0. One can speculate that the appearance of the 
global action of a compact group is related to the presence of a "generic" singularity 
but there is no general proof yet. Note that the generic Neumann system has no 
symmetries given by a compact group. 

The singularities of the spectral curve appeared in two different roles in our study. 
The singularity that is a consequence of the confluency is "generic" in that it appears 
uniformly for all values of the energy-momentum map. The rotational symmetry 
shows as the extension by complexified group of rotations C* defining the generalized 
Jacobian that appears globally. The "sporadic" singularities, which correspond to 
the singular values of the energy momentum map (see proposition I3.8fl are strictly a 
local phenomenon. In those cases the extension by C* and the resulting rotational 
symmetry does not extend globally. Algebraically speaking both singularities are the 
same, but the "generic" singularity appears globally and thus give rise to a rotational 
symmetry. Sporadic singularities on the other hand appear when the level sets of 
energy momentum map are singular (orbits of lower dimension, heteroclinic and 
homoclinic orbits). It would be interesting to describe the isospectral sets of the 
singular spectral curves. Note that when we introduced generic singularity we made 
sure that we used the subset of the singular isospectral set, consisting of regular Lax 
matrices. In the study of sporadic singularities, non regular part of the isospectral 
set should not be avoided. It is our conjecture that the singular isospectral sets that 
induce homoclinic or heteroclinic orbits should pose an obstruction to the existence 
of global action of compact groups. 

In a somewhat more ambitious and speculative vein, one could study the 
relationship between symmetries of certain PDE's and generic singularities in 
appropriate spectral curves of infinite genus. Maxwell-Bloch equation for example 
can be viewed as a chain of confluent Neumann systems [18j [19], whose symmetries 
indeed reflect in a symmetry of the whole Maxwell-Bloch system [20]. The description 
of Maxwell-Bloch system with generalized Jacobians of singular spectral curves of 
infinite genus should be to some extent analogous to our results. 
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